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Abstract. The main goal of this work is to prove that every non-negative 
strong solution u(x,t) to the problem 

u t + {-A) a/2 u = for (x,t) 6 R n x (0,T), < a < 2, 
can be written as 



u{x,t)= P t (x -y)u(y,0)dy, 



where 

P t (a;) = — Lp 
and 

P(x) := f e ix <-^"d^. 
JtL" 

This result shows uniqueness in the setting of non-negative solutions and ex- 
tends some classical results for the heat equation by D. V. Widder in 1 1 21 to 
the nonlocal diffusion framework. 



1. Introduction 

The heat equation has the hyperplane t = as a characteristic surface and this 
causes that the Cauchy problem with initial data on t = is not well possed in 
general. However there is a clear agreement between the Principles of Thermo- 
dynamics and the model of transfer of heat given by such an equation. This is 
reflected into the fact that the initial temperature evolves in time as the convolu- 
tion with a kernel, giving rise to an average, and smoothing out in this way the 
potential effect of sharp thermal differences (this is in agreement with the entropy 
effect of the Second Principle of Thermodynamics) ; also uniqueness holds true un- 
der a positivity assumption on the function (this is in agreement with the so called 
Third Principle of Thermodynamics, according to which temperatures are always 
positive if measured in the Kelvin scale). In this sense D. V. Widder in [12] proved 
the following classical result: 

Assume that u : R n x [0, T) C Tf ->• M is so that 

u(x, t) > 0, u € <tf(R n x [0, T)), m, u x%Xz E %f(R n x (0, T)), 
and satisfies 

u t (x, t) - Au(x, t) = 0, (x, t) e R™ x (0, T), 
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in the classical sense. 
Then 

l f ( — \x — y\ \ 

u ( x > t) = ta / u (y> °) exp ( — t* — dy - 

(47Tt)2 J Rn \ to J 

In this paper we obtain a similar result for the nonlocal heat equation, that is 
the equation, 

(1.1) u t + (~A) a/2 u = for (x, t) e R n x (0, T), < a < 2, 

in which the diffusion is given by a power of the Laplacian (we refer to [5] and [TT] 
for basic definitions and properties of the fractional Laplace operator). 

It is worthy to point out that for every aG (0,2] the operator in (11.11) does not 
satisfies the so called Hadamard condition, that is, the Cauchy problem is ill possed 
(see [7] for more details). As a consequence we face here a problem similar to the 
one of the heat equation. 

It is well known that the operator in (jl.ll) and its inverse are nonlocal. In a sense, 
equation (11.11) is perhaps one of the simplest classical examples of pseudodiffcrential 
operators. Notice that if we consider the problem 

u t + (~A) a/2 u = for (x, t) E R n x (0, T), < a < 2 

with initial datum u(x,Q) = uq(x) assuming, say, uo <= ^(K™) H L oa (R n ), then a 
solution is obtained by the convolution with the kernel 

(1-2) Pt(x) = -n-P(-vr) , 

where 

(1.3) P(x) := f e fa *-l* |0 d£. 
That is, a solution is given by 

(1.4) u(x,t)= P t (x-y)u (y)dy. 

We look here for a class of solutions of the fractional parabolic equation, such that 
the sign condition u(x, t) imply that u, necessarily, is of the form (|1.4p with 
Uq(x) replaced by the trace u{x, 0). This is the type of extension that we propose 
for the classical result of Widder to the nonlocal equation in (|1.1[) . 

In this article we will describe several interpretations of what solution to equation 
(jl.ip means, and consider weak, viscosity and strong solutions, in a sense that we 
now make more precise. 

Weak solutions. Consider the space 

if Q / 2 (M") := [u : W -> R measurable : [ \ u f> x )\ dx < oo 
endowed with the norm 

ll u l|j^°/ 2 (R") = 



\u(x)\ 
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If u G Jzf Q / 2 (R") then (-A) Q / 2 u can be defined in the weak sense as a distribution, 
that is we can compute the duality product ((— A) Q / 2 it, tp), for every ip in the 
Schwartz class. 

Definition 1.1. We say thatu(x,t) is a weak solution of the fractional heat problem 

u t + (-A) Q / 2 u = for (x, t) G E" x (0, T), 
u(x, 0) = uq(x) mR n , 

if the following conditions hold: 

i) u G L 1 ( [0, T'] , Jz? a / 2 (R n ) ) /or every T < T. 

ii) ttefiT),^!")). 

iii) For every test function ip G Co°(R" x [0, T)) and < T" < T one /ias i/iai 



(1.5) / / -u(x,t)^ t (x,t) + u(x,t){~A) a/2 ^(x,t) 



dx dt = 



u(x,T r )ip(x,T') dx — I uo(x)>p(x, 0) dx. 

Condition ii) is imposed so that the right hand side of equality (|1.5[) makes sense 
for every T'. Notice that the continuity would follow in any case from the left hand 
side of (|1.5[) due to the integrability condition in i). 



Viscosity solutions. Consider Q = R™ x (0,T) and 

K' 2 (Q) = {/ : Q -> R I /t G ff(Q), fxi,xj G tf(Q) 

and sup |/0M)l < C(l + |x|) p } 
te(o,T) J 

Definition 1.2. ^4 function u G ^(Q) is a viscosity subsolution (resp. supersolu- 
tion) o/ 

(1.6) Ut + (-A) Q / 2 ( u ) =0 

in Q if for all (x 7 t) G Q and tp G (Q) such that u — tp attains a local maximum 
(minimum) at (x,t) one has 

p t (x,i) + (-A) a / 2 (ip(x,i))) ^0 (resp. ». 

M^e say f/ia£ zt G c €(ff) is a viscosity solution of (|1.6[) m Q is both a viscosity 
subsolution and super solution. See [6] for more details about these type of solutions. 



Strong solutions. We know (see Proposition 2.1.4 of [TDJ) that if u G Jz? a/2 (R n )n 
^l„^ 7 (M n ) (or if a > l) ; for some 7 > then (-A) Q / 2 u is a continuous 

function and it is well defined as the following principal value 

(-Ar^u(x) := C(n, a)P.V. / u{x) - "ff dy, a G (0, 2). 

jb" f — y\ 

Here C(n, a) denotes the constant satisfying the identity 

(-A) Q/2 u = ^- 1 (\^\ a ^v), (el",«ey,«G (0, 2); 
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that is, 

C(n,a)=([ i^^r 1 



see 0. 



In order to allow a larger class of functions, and following standard procedures 
in the theory of singular integrals, we will define hereafter the principal value as 
the two-sided limit 

f u(x)-u(y) I u(x) -u(y) 

V \x - y\ n +« * e^o J \ x - 

{y I e<\x-y\<l/e] 

We will use the same notation for this extended operator. When u € «5f"/ 2 (R™) 
this definition coincides with the usual one. 

Definition 1.3. We say that u(x,t) is a strong solution of the fractional heat 
equation 

u t + (-A)"/ 2 u = for (x, ()el"x (0, T), 
if the following conditions hold: 

i) u t (x,t) £ %f(R" x (0,T)). 

ii) u e "af(R" x [0,T)). 

iii) The equation is satisfied pointwise for every (x,t) £ R™ x (0,T), that is, 

u t (x, t) + C(n,a) P .V. [ ^ ^ dy = 0. 
Note that if u(x,t) is a strong solution then 

Observe also that if u(x,t) is a strong solution of the fractional heat equation 
satisfying u(x,t) G L 1 ([0, T'], Jz? Q / 2 (R™)) for every T < T, then u(x,t) is a weak 
solution of the fractional heat equation (as a byproduct of our results, we will see 
that this holds true for non-negative strong solutions, see the forthcoming Corol- 
lary [SU). 



Our main contribution in this setting is the following Widder's type result 
Theorem 1.4. If u is a strong solution of 

u t + {-A) a/2 u = for 0, t) 6 R" x (0, T), < a < 2, 

u(x,t) = P t {x- y)u(y, 0) dy. 



Note that, if for a given polynomial p with real coefficients we consider the 
differential operator p(D), then we obtain that 



p{D)P t (x)= / e^ptiOe-*"'! d£- 

Since e - *^" is a tempered distribution we deduce that G C°°(R" x (0, oo)) (see 
for instance [5] and [5] for more details). In particular, and as a consequence of 
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Theorem II. A\ we get that if u is a non negative strong solution of the fractional 
heat equation then u € C°°(K™ x (0,T)). 

The paper is organized as follows. In Section [5] we prove a uniqueness result 
for weak solutions that, in turn, will be the key step to obtain our representation 
theorem. Next in Section [3] we prove the main result for strong solutions. We start 
by proving a comparison result that will allow us to show that every positive strong 
solution u{x, t) is bigger than or equal to the convolution of the trace u(x,0) with 
the kernel Pt(x). By a scaling argument we will prove that any positive strong 
solution is also a weak solution and then, by the uniqueness result of the previous 
section, we will conclude with the proof of the theorem. 

Finally, in the last section we establish the pointwise behavior of positive strong 
solution that has an interest on its own and that could provide and alternative 
proof to our representation result. 

2. Uniqueness for weak solutions 

To start with, we prove a uniqueness result for weak solutions with vanishing 
initial condition. The proof is quite complicated and it involves many fine integral 
estimates. The nonlocal feature of the problem also makes localization and cutoff 
arguments much harder than in the classical case. 

Theorem 2.1. Set T > 0, a £ (0, 2) and let u be a weak solution of the fractional 
heat equation 

u t + {-&) a/2 u = for (x, t) e M™ x (0, T), 
u(x,0)=0 inW n . 

Then u(x,t) = for every t G (0,T) and a.e. x £ M. n . 



(2.1) 



Proof. We must show that u(x,to) — for an arbitrary to £ (0, T) and x S R". 
For this, we fix R > and 9 € ^^{B^) and we will prove that 

/ u(x, to)O(x) dx = 0. 
it" 

For any t £ [0,to), we define 

<p(x,t) := (0(-)*Pt o -t(O)(z), 
where Pt is the kernel defined in (|1.2j) and (|1.3[) . By [5] (see also [5] and [1]) we 
know that 

11 C 
( 2 - 2 ) c 1 + \x\ n + a ^ P{X) ^ l + \x\ n +° 

Therefore 

Since 

we have that 
(2.3) 

Now we claim that 

(2-4) \<p{x,t)\ = \9{x)*P^ t {x)\ < 1+ ?\ n+a -. 



= ierw,*), 

ip t - (~A) a / 2 ip = for (x, t) € M" x [0, to), 
(p(x,to) = 9(x) inM™. 
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where C\ depends of n, a, Rq, M := ||^||zi«>(b r ) and to. 

Indeed, considering without lost of generality that 2R > 1, we will distinguish two 
cases: 

When \x\ ^ 2i?o, using that Pt is a summability kernel in L 1 we have that 



p t o -t(y)0(x - y) dy 



(2.5) 



s$ M 



< M 



P t0 - t (y)dy = M 



1 + (2R ) n+a 
1 + \x\ n+a 

Cl 



1 + 1x1"+" 



where c\ = Ci(n, a, Rq, M). 

Consider now \x\ > 2Rq. Note that from (11.21) and (12.21) . we obtain that 



(2.6) 



P t (y) «S 



C 



f- 



Ct 



\y 



Then, since for \x — y\ ^ R it follows that \y\ ^ we have 

2 n+a C(t Q - t) 



Pt -t(y) < 

As a consequence, 

P to -t(y)0(x - y) dy 



\n-\-a. 



Pt o -t{y)0{x - y) dy 



x-y\^R 



(2.7) 



< 2 n+a CM\B Ro 

< 2" +a 2CAf \B E 



*0 



1 + \x\ n+a 



C-2 



1 + \x\ n+a ' 

where C2 = 02(11, a, i?o, Af, to) and |-Br | denotes, as usual, the Lebesgue measure 



of the ball. Hence, (j2l| follows from ([23)) and ([2T7]) . 

Also, applying ()2.7j) to the derivatives of G ^"(-Br,,), we also have 

C 2 



(2.8) |V^(x,<)| = \V6(x) * P t0 -t(x)\ < r+ +n . 

where C2 depends also of n, a, i?o, M' :~ \\ \J6\\l^(b Rq ) and to- 
Then, from ([231) and the fact that u € ^([O, i ], Jz? q/2 (R")) we deduce that 



(2.9) 

Let now <j) e c if° 
(2.10) 



M 2 := 



In 



\u(x, t)ip(x, t)\ dx dt < 00. 
be such that 

XB 1/2 < <j> < XBi- 
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For R > 2R we define 4> R (x) = (f> (-|) and 

ip(x,t) := (p(x,t)(p R (x), 

as a test function of problem (|2.1|) . Also, by an explicit computation, one sees that, 
for any f,ge Jz? Q / 2 (R n ), 

(-£) a ' 2 (fg)(x) = f(x){-/\r' 2 g{x)+g{x){~^rl 2 f{x) - B(f,g)(x) 
with B(f,g) the bilinear form given by 

(f(x)-f(y))(g(x)-g(y)) 



B{f,g)(x) :=C(n,a) 



\x - y\ n+a 



dy. 



Appling this formula (for a fixed t) to the functions tp and 4> R and recalling (|2.3p 
we obtain 



(2-11) 



(-A)"/ 2 ^ - (^(-A) Q / 2 fl + fan - B(cp, <t> R ). 



Moreover, the function uip evaluated at to is u(x, to)6(x)(j> R {x), thanks to the ter- 
minal time condition in (j2.3[) . Thus, considering ifi as a test function in (|2.1[) and 
using that ip is a solution of the problem (j2.3[) we obtain 



(2.12) 



u(x, to)0(x)4>n(x) dx 



dx dt 



uip t (x,t)faix) - u(-A) a / 2 (^{x,t)) 
uB(cp, (j> R )(x, t) - utp(x, t)(-A) a / 2 <j> R (x) 
wp(x,t)\ \ (-A) a/2 <t> R (x)\dx dt 

\u(x,t)\ \B((p,<j) R )(x,t)\dxdt 



dx dt 



Then, since 9 is supported in B Ro and Rq < R/2, and recalling (|2.10[) . we conclude 
that 

u(x,t )9(x)dx < / / \u(p(x,t)\ \(-A) a/2 cf) R (x)\dxdt 



(2.13) 

It remains to show that 



\u(x,t)\ \B(tp,(j> R )(x,t)\ dxdt 
h(R) + C(n,a)I 2 (R). 



lim h(R) + I 2 (R) = 0. 

R— >oo 



Indeed, since 

K-A) Q / 2 Mz)i=^ Q |((-A) Q/ y) ( 

by (|23)1 . it follows that 

7i(i?) «C C ir a / ° / \uip(x,t)\dx 
Jo JR" 



^ Cni? Q , 



dt < C M 2 R~ a 
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and so 
(2.14) 



lim h(R) = 0. 



Now we are going to estimate I2(R)- For this we cover R" x l n with six domains 
suitably described by the radii -R/4, R/2, R and 2R and represented (for n — 1) in 
the following picture: 



2R 




Therefore 



r,2/i 



where 

A 1 := {(x,y) : |x| > R/2, \y\ ^ R/4}, A 2 := {(x,y) : \x\ < i?/4, | v | > R/2}, 
A 3 := {{x,y) : \x\ > 2E, < |y| < fl}, A 4 := {(x,y) : J2/4 < |z| < R, \y\ > 2R}, 
A 5 := {{x,y) : R/4 < \x\ < 2R, R/4 < \y\ < 2R} 



and 



C := {(x, y) : |x| < R/2, \y\ < R/2} U {(x, y) : |x| > R, \y\ > R}. 
From ()2.10j) . we know that </)r(x) — <t>R.{y) — if (x, y) 6 C, and so 



rt 



I 2 (R) = / / \u(x,t)\ 
Jo Jr™ 

(2.15) < E^W, 

where 



k=l 



\ip(x,t) - ip{y,t) \ \cj> R {x) - (j> R {y)\ 
\x - y\ n+a 



dy dx dt 



li"(R) = 



for k = 1, . . . , 5. 



J At 



, / .si \f{x,t) - ip{y,t) \ \<j> R {x) - <j) R {y)\ 

\u(x, t) | _ ^ n+a dy dx dt, 
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We are going to estimate each of these five integral separately. For (x, y) G A\ 
we get that \x — y| ^ C\x\. Moreover by ()2.4j) it follows that 

|y(M)l + lp(g,*)|< 1 + |y|n+a - 

Therefore 

(2.16) < cf /" M^l^dt. 
Following the same ideas, since for (x, y) G Ai we obtain that 

(2.17) b(^t)l + lv(y»*)l<Tq^» and l^-y|>c|y|, 



then 



JO ik|^fl/4 Fl J\y\>R/2 Wl 



C 



(2.18) < CUT"/ I T^dxdt. 



Jl! r l u(x,t)\ 

J\x\^R/4 



Also, since for (x, y) G A3 (|2 . 1 7[) is satisfied, then 



(2.19) 



If* {R) SCR- f° I ^$d xd t, 

JO J\x\^2R \ x \ 



Similarly, using again the good decay of if and the fact that \x — y\ ^ C\y\ for every 
(cc,y) G A.4, we obtain that 



(2.20) lti R ) < CR ~ a I ° I 

Jo J h 



Hx ^dxdt. 



R/A<\x\<R Fl 

Then, using the Monotone Convergence Theorem, by (|2.9I) and the fact that u G 
L 1 ([0,t o },^' a/2 (M. n )), from (l2~16l) . (l2~18l) . (l2~19l) and (l2~20l) it follows that 

(2.21) lim l£ x {R) = lim I^iR) = hm I 2 A3 (i?) = lim I?*(R) = 0. 

R— >oo R— >oc _R— >oo R— >oo 

To estimate I^ 5 (i?) we will treat separately the cases a G [0, 1) and a G [1, 2). 

We start with the case a G [0, 1) and (a;, y) & A 5 . Since in A 5 the roles of x and 
y are symmetric, we deduce from (|2.4|) that, in this case, 

(2-22) |^(i,t)| + |^,t)|<-— . 

Also 

C 

|<Mz) - <My)| < -^k - y|. 

Thus 

(2.23) I^(R)^Cf tn f J!MI / 
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By the change of variables y := x — y, it follows from (|2.23|) that 



RJo Jjl^^R \x\ n+a jK^\ x -y\^2R\y\ 

< — [ to ( KM)! f 1 



<|x|<2fl l x l ^|£|<4fl 

|u(a;, i)| 



(2.24) < CiT*/ / L-S^didt. 



Therefore, using that u G L 1 ([0,to] ) -^ a/2 (K n )) ) we conclude that 



(2.25) lim l£ 5 (R) = 0, when a E [0, 1). 

R— >oo 



We consider now the case a £ [1,2). By (|2.8[) . we get that 



(2-26) < i + uin+> -^ 



for some z in the segment joining x and y. We take the family 

R 



Q := <j (i,j)£A 5 : \x-y\ < 

Note that, if (x, y) £ Q then every point 2 lying on the segment from a; to y verifies 
\z\ > C|x|. Hence, (|2.26p and the previous estimate for (j>n, gives that 



A) J(x,y)£Q \X—y\ n 



(■lit 



* Jo l, 9) ^ (l ' i) '«n^ ¥ldt 

C /" to /" |u(ar, £)| f 1 , , , 

< ™ / / / ... u^ dydxdt 



RJo J^\x\^2R \A n+a J^\y\^R\x-y\ n+a - 2 

c f to r \u(x,t)\ r i 
*" W J^\x\&r \*\ n+a J%<\x-m™ \v\ n+a - 2 y x 1 

(2.27) < OR'-* f f l T^dxdt. 

Jo J&^\x\42R \ X \ + 

On the other hand, if (x, y) E A$\Q we have that 



(2.28) 
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Then by (|2~22]l and (f2~2"8]) it follows that 



\(4>r{x) - (j> R (y))(ip (x,t) - <p(y,t))\ 
J(x,y)eA 5 \Q \ x - Ij\ 

C f t0 f , , ^.\{<p(x,t) - ip(y,t)) 



u(x, t) | j \ n +a ' ^"1/ dx dt 



, , \u(x,t)\ ; - : : \ dv dx dt 

u JO J(x,y)eA 5 \Q \ x V\ 



C f to [ \u(x,t)\ 1 , , , 

sC — / / , \ ' n — — - — T dydxdt 

RJo J( x .v)eA,\Q \x\ n+a y ^- 1 



RJo J(x,y)eA 5 \Q \A n+a \VY 

< c f t0 f \u(x,t)\ r i 

(2.29) < CR- a f a [ \^S.dxdt. 
Therefore, from (|2T2Tf and (g^]) 



to 



- <p R (y))(<p(x,t) - <p(y,t))\ 



lt(R) < / / | U ( g ,t)| I^vv rfj^m ™,-,n dydxdt 

JO J(x,y)GQ \ x — V\ 



to 



(2.8(1. + CiT" f° i l -^^dxdt. 



(x) - 4> R {y))((p(x, t) - ip(y, t))\ 



J§^\x\^2R 



Then, since u 6 i 1 ([0, to], Jz?"/ 2 (R™)), using the Monotone Convergence Theorem 
we obtain 

(2.31) lim I^'{R) = 0, when a £ [1,2). 
That is, by ([235]) and (j2~3T]) . we get 

(2.32) lim |/^ 5 (i?)| = 0, whenever a G (0,2). 

R— >oo 

Putting together (|2~T51) . (pT2l) and (f232l it follows that 

(2.33) lim I 2 (R) = 0, when a e (0,2). 

R— >oo 

Therefore, from (f2"33]l . by flUTD and ([233]) we conclude that 



lim / u(x, to)9(x) dx = 0, 
for an arbitrary 9 £ ^(B Ro ), R < 2R. □ 



3. Uniqueness for strong positive solutions 



In this section we will establish the representation of the positive strong solutions 
of the fractional heat equation as the Poisson integral of the initial value. That is 
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Theorem 3.1. Let (x,t) £ K™ x (0,T). If u{x,t) ^ is a strong solution of the 
fractional heat equation, then 

P t (x - y)u(y, 0) dy = u(x,t). 

To prove this theorem we will need some previous results that we present as 
follows. First of all, we establish that, among all possible positive solutions of the 
fractional heat equation, the smallest one is given by a formula that involves the 
convolution with the fractional heat kernel (see Lemma 13.31) . To prove it we will 
use the following 

Lemma 3.2 (A maximum principle). Set Dt := O X (0,T) and let v(x,t) £ 
^(Q x [0,T)) satisfy, pointwise, the following problem 



(3.1) 



v t + (~A) a / 2 v ^ for (x, t) £ D Tl 
v(x,t)^0 in (R n x [0,T)) \D T . 



Then v < in Q X [0,T). 

Proof. Fixing an arbitrary T' £ (0,T), we define 

v{xq, to) '■= _max v(x,t). 

Ox [0,7"] 

Our goal is to show that v(xo,to) ^ 0. The proof is by contradiction, assuming 
that 

(3.2) v(x ,t )>0. 

Then, (xo,t ) cannot lie in (dfl x [0,T)) U (f2 x {0}), since v ^ there, thanks to 
the boundary conditions in (|3.1|) . As a consequence, (xo,to) nes in ^ x (0, T'] and 
therefore v t {xQ,to) = 0. Therefore the equation in (|3.ip implies that 

/ v(*o,to)-r>(l/,to) 
' 7 H n\n \y-x \ n + a 

Since v(y,to) < for y S R n \£l, thanks to (|3.1[) , we obtain that the latter integrand 
is strictly positive, due to (|3.2[) . and this is a contradiction. □ 

Now we are able to prove that strong solutions are upper bounds for the kernel 
convolutions: 

Lemma 3.3. Let (x,t) £ M" x (0,T). If u(x,t) ^ is a strong solution of the 
fractional heat equation then 



(3.3) /:= / P t (x-y)u(y,0)dy^u(x,t), 

JK™ 

where Pt(x) is the function defined in (|1 .2[) . 
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Proof. First of all, observe that the integral I = I(x,t) exists, for is given by the 
integration of two (measurable) positive functions, although we do not know a priori 
that I is finite. However, this will be a consequence of our result that gives the 
inequality / ^ u(x,t) for every (x,t) G R™ x [0,T). To this aim, we let 

r i, izKiz-i, 

(3.4) <t> R (x):=l R-\x\, R-l^\x\^R, 

[ 0, \x\ > R. 

We define 

v R (x,t):= [ P t {x-y)<t> R {y)u{y,0)dy = {P t {-)*<t>Ru{;Q)){x). 

Then 

+ {-A) c '/ 2 v R = for (x, t) e E" x (0, T), 

V R (x, t) > for (x, t) e M" x (0, T), 

VR (x, 0) = <f> R (x)u(x, 0) in R™. 

Let \x\ > R. As u(x,t) S ^(R" x [0,T)) we can define the real number 

M R := sup u(y,0) < oo. 

By (|2.6p we have that 

< v R (x, t) < M R [ P t {x- y) dy 
Jb r 

< CM R [ j—^dy 



^ C(T,M R ) f 

JE 



= C(T,M fl ,n) 



\x\ - R\ n + a 

for any (x,t) G (R™ \ Br) x (0,T). 
Then, for every e > it follows that 

(3.5) 0^v fl (:M)^£, for any \x\ > p, i € (0,T), 
where 

Moreover, as u(x,t) > in R™ x [0,T) we obtain that 

(3.6) < v R (x,t) < e s$ e + u(x,t), for any |x| ^ p,t E (0,T) 
and 

(3.7) ur(x, 0) = 4> R {x)u{x, 0) ^ w(a;, 0) ^ e + 0), for any \x\ ^ p. 
Consider the cylinder 

D p ,t = B p x (0,T). 

We define the function 

w(x, t) := v R (x, t) — u(x, t) — e. 
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Then, by flU) and IpTf]) . we get that w(x,t) ^ in R™ x [0,T) \ L> P , T . Therefore, 
since u(x,t) is a strong solution of the fractional heat equation in R™ x [0,T), 
applying Lemma 13.21 in -D p ,t to the function w{x,t) we have that 

v R (x, t) ^ e + i), for |x| ^ p and i G [0, T). 
Therefore, from (|3.6j) . it follows that 

u fl (x,t) < e + u(x,t), for x £ R" and i G [0,T). 
Since e is fixed but arbitrary, the previous inequality implies that 

vn(x,t) < u(x,t), for every x G R" and t G [0,T). 

Finally, by the Monotone Convergence Theorem, as lim <fin = 1, we conclude that 

R— >QO 

0^v(x,t) = lim VR(x,t) = / P t (x - y)u(y,0)dy ^ u(x,t). □ 

R^OC Jm,n 



By a simple time translation, we obtain from Lemma 13.31 the following: 

Corollary 3.4. Let < r < T and (x,t) G R" x (0,T - r). If u(x,t) ^ is a 
strong solution of the fractional heat equation then 

(3.8) J P t {x - y)u(y 7 r)dy ^ u(x,t + t), 

where Pt(x) is the function defined in (|1.2|) . ^4s a consequence, for every x G R™ 
and t 6 [0, T — r) we ftawe 

(3.9) / / P t (x-y)u(y,T)dydT^ / w(x,i + r)dr. 



Moreover we have the following 

Corollary 3.5. Let (x,t) G R" x (0,T). If u(x,t) ^ is a strong solution of the 
fractional heat equation then, u(-,t) G Jz?"/ 2 (R™). 

Proof. Let < T < T. Taking t = T — T'm, §n$, from (J22|), we get that 
T-T' f u{y,T) 



C V (T - T') 



\x-y 



dy < / Pt~T'{x - y)u{y,r)dy 



(3.10) < u(x,T-T' + t) < co. 

Let 

1 1 



Then, since 



C(T) := min < 1, - 

1 (T - T')^r 



(T-T") 2 ^ 2 + \y\ n+a 'l + \y\ n+a ' 

taking x = in (I3.10p . we conclude that 



/ t; 1 i . n dy < oo. D 

V (1 + M n+Q ) 
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We also have 

Corollary 3.6. Let (x,t) G M" x (0,T). If u(x,t) ^ is a strong solution of the 
fractional heat equation, then u(-,t) G L 1 ([0,T'],^f a / 2 (R n )) for every < T' < T. 

Proof. Take an arbitrary < V < T. By (EU with i = T - V we have that 

u(x. t)6lt. 



\ \ P T - T >(x-y)u(y,T)dydT < / 

Jo JR™ JO 



Then, as u G ^(ff 1 x [0,T']) we get 

fT' 



/ / P T -t' (x - y)u(y,r)dydT < oo. 

JO JR™ 

Therefore doing as the proof of Corollary 13.51 we get that 



cty <m < oo. 



(1 + | y |n+a) 

That is, u G ^([O, T'], if Q / 2 (R n )) for every < T < T. □ 

Note that Corollary 13.61 affirms that if u(x, t) > is a strong solution of the 
fractional heat equation then u is also a weak solution of the same equation. 

Now we are able to prove our main result: 

Proof of Theorem\Tj\ By Corollary EH we get that u G L 1 ([0, T'\, Jz? Q / 2 (R")) for 
every < T' < T. Moreover if we define 



p(x,t) := P t (x-y)u(y,0)dy, 

JR™ 

by (EU) with t = T-T',we also have that p G L x ([0, T'], «^ Q / 2 (R")). Let now the 
function 

w(x,t) :— u(x,t) — p(x,t) 0. 
It is clear that w is a strong solution of the fractional heat equation. Moreover, as 
w G i 1 ([0, T'], Jz? a / 2 (R n )), then w(x,t) is also a solution in the weak sense with 
zero initial datum. Therefore applying Theorem 12.11 we conclude that w(x,t) = 
for every (x, t) G R" x [0, T). □ 

3.1. About viscosity solutions. As was said at the beginning of this work, it is 
natural to consider viscosity solutions of the fractional heat equation. Our purpose 
here is to describe some cases in which a positive viscosity solution has the unique 
representation in terms of the kernel Pt ■ 

Proposition 3.7. Let {«„} be a sequence of non-negative, strong solutions of the 
fractional heat equation converging, uniformly over compact sets, to a given function 
u. Then u is a viscosity solution of (II. ip satisfying 

(3.11) J P t (x-y)u(y,0)dy^u(x,t), (x, t) G R" x (0, T) 
and 

(3.12) !ieI 1 ([0,T'],if a / 2 (l n )) for every < T' <T. 
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That is, the conclusions of Lemma \3.3\ and Corollary 1 3. 6\ are satisfied. 
Proof. By Lemma [3~3l it follows that 

P t (x - y)u n (y, 0) dy < u n {x, t), (x, t) e R" x (0, T). 



Applying the Fatou Lemma we obtain (|3.11l) . Therefore, doing as in the proof of 
the Corollary 13.61 we conclude (|3 . 1 2[) . Note also that, by the comparison principle 
(Corolary 2.1.6 of jTD]), u n ^ is a viscosity solution of (|1.1[) for every n e N. 
Therefore, since u is the uniform limit over compact sets of viscosity solutions, we 
get that u is also a viscosity solution of (jl.ll) . □ 

Notice that to conclude the equality in (|3 . 11 1) we would need to know that u is 
a weak solution. 

If we add a monotonicity condition over the sequence {u n } then we obtain the 
following 

Proposition 3.8. Let {u n } be a monotone sequence of non-negative, strong solu- 
tions of the fractional heat equation converging, uniformly over compact sets, to a 
given function u. Then u is a strong solution of (|1.1[) satisfying 

(3.13) / P t (x-y)u{y,0)dy = u(x,t), (x, t) e R" x (0, T). 

JR™ 

Proof. Since, for every n £ N, % ) satisfies Theorem \1A\ by the Monotone 
Convergence Theorem we obtain (|3.13[) . Clearly this implies that u ^ is a strong 
solution of qi.ip . □ 

Remark 3.9. It would be interesting to find the biggest class of positive viscosity 
solution for which the representation property (|3.13|) holds. This seems to be an 
open problem as far as we know. 

4. Further results. 
As in the local case (see [13]), given a solution u, one can define the enthalpy 

v(x, t) = u(x, s) ds 
Jo 

which is also a solution of the fractional heat equation. Indeed we have the following 
result: 

Lemma 4.1. Let [x, t) G R" x [0, T). If u(x, t) is a strong solution of the fractional 
heat equation with vanishing initial condition, then the entalpy term 

v(x, t) — u(x, s) ds 
Jo 

is also a strong solution of the fractional heat equation. Moreover if u is positive 
the function v is increasing in t for x fixed and a/2-subharmonic as a function of 
x in the weak sense. 

Proof. Let (x, t) € R™ x [0, T). First of all note that v(x, t) satisfies the conditions i)- 
iii) of the Dcfinition ll.3l Therefore, since u(x, t) is a strong solution of the fractional 
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heat equation, by the Fundamental Theorem of Calculus and Fubini Theorem it 
follows that 

C(n,a)P.V.[ VM dy 

JR™ \ x ~ V\ 

= C(n,a P.V. / -y : dy 

J«» f - y\ 

= C(n,a P.V. / / — : — — — dsdy 

Jmn Jo \ x - V\ + 

(-A) a/2 u(x,s)ds 
i 



u s (s, t) ds 

o 

= — u(x, t) 
= -v t (x,t), 

for any (x, t) G R" X (0, T]. Then, v(x,t) satisfies the fractional heat equation in 
the strong sense. Also if u ^ from the calculations done before we deduce that 

-(-A) a/2 v(x,t) = v t {x,t) = u{x,t) Ss 0. 

So v is a/2-subharmonic as a function of x and increasing in t for x fixed. □ 

Lemma |4~T1 shows that the enthalpy belongs to the special class of positive strong 
solutions that are also a/2-subharmonic. This class naturally satisfies a polynomial 
estimate, as shown by the following result. 

Proposition 4.2. Let (x,t) el"x [0,T). Let u(x,t) ^ such that 

i) u(x,t) is an a /2-subharmonic function with respect to the variable x. 

ii) u(x,t) is a strong solution of the fractional heat equation. 

Thenu(x,t) C(t)(l + \x\ n+a ) if (x,t) G R™ x [0,T). 

Proof. By hypothesis it is clear that u is an a/2-subharmonic function with respect 

to the variable x for t fixed and therefore u is increasing in time. 

Let now < t 1 < T and < t < T — 1\. By Corollary UO] we have that 



Pt(x — y)u(y, t\)dy ^ u{x, t + for any < t < T — 1\. 
Therefore 

M to := / P to (y)u(y, ti)dy < u(0, t + h) < oo. 
Our objective is to show that 

(4.1) \u(x,t 1 )\^C(l + \x\ n+a ). 

Once this is done, using that u is increasing in time, we would get 

s$ u{x,t) sC u(x,h) ^ (7(1 + \x\ n+a ) for every (x,t) 6 R n x (0,ti). 
But, since io and £i are fixed but arbitrary, we would conclude that 

\u(x,t)\^C(l + \x\ n+a ). 
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So, we are left to showing that (|4.1[) is true. Note that from Corollary 13.51 we have 
that u S J?f Q / 2 (]R"). Then as u is a/2-subharmonic, by Proposition 2.2.6 of [10] u 
also satisfies the following mean value property: 
(4.2) 

u(xo,t) ^ / 7a (y — x a )u(y,t)dy for every xq e £1 C E n and A ^ dist(x , <9f2), 
where 

lx (x) := (-A) a/2 r A (x), 

and 

Here T is a C 1 ' 1 function that coincides with $(x) := c|a;| a n outside the ball B\ 
and r is a paraboloid inside this ball. Note that $ is the fundamental solution of 
(_A)"/ 2 . From flU) we have 

u(z I fi)(l + M" +c T 1 ^ {l + \x\ n+a )- x f 7A(l/)«(a:-l/,*i)di/ 

= (l + |xr+«)- 1 / lx (y)u(x-yM)dy 



'{y-\y\»-} 

+ (1 + / 7x(»)«(»-y,*i)dl/ 

(4.3) := 7 1 (a;)+/ 2 (x). 
Choosing 

(4.4) A=|x|/4, 

we have that \x — y\ ^ 5|y|. Therefore, by Proposition 2.2.3 of (TU], we obtain 

u(x - y,t x ) 

'{w-\y\>>>} 

u(x - y,h) 
\x - y\ n+a V 

(4.5) < C\\u\y a/HRn) -.= 0,. 

Moreover, using that the fractional Laplacian of a paraboloid is bounded, by (12.2 
and (14.41). it follows that 



^ c(i + i^r+")- 1 



h{x) ^ C{\ + \x\ n+a )~ 1 / u{x-y,t x )dy 

J{y.\y\<\} 

< c(i + M"^)- 1 / u(z,h)dz 

J{z: \x-z\^\} 

^ C{l + \x\ n+a )- 1 j u {zM)Pi{z)-^—dz 

J{z: |*|<2|x|} "l( 2 J 

< C(l + |a;|" + ")- 1 (l + (2|x|)" +Q ) / ufotO^i 

J{z: |z|^2|a:|} 

< CMi(l + |a;|»+«)-i(i + {2\x\) n+a ) 
(4.6) < C(n,a)u{0,l + h):=C 2 . 

By l|4.3|) . (|4~5l) and (l4~6l) . we obtain (|4~T|) and we conclude the proof. □ 
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Remark 4.3. It would be interesting to prove that the pointwise behavior in Propo- 
sition \4~J^ and some comparison arguments provide an alternative proof to our main 
result, namely, the representation of every solution as a convolution with the asso- 
ciated Poisson kernel. 
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